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Q\ [ We study two-dimensional Dirac fermions in a random non-Abelian vector potential 
by using lattice regularization. We consider U(N) random vector potential for large 
N. The ensemble average with respect to random vector potential is taken by us- 
ing lattice supersymmetry which we introduced before in order to investigate phase 



structure of supersymmetric gauge theory. We show that a phase transition occurs 
at a certain critical disorder strength. The ground state and low-energy excitations 
are studied in detail in the strong-disorder phase. Correlation function of the fermion 
local density of states decays algebraically at the band center because of a quasi-long- 
range order of chiral symmetry and the chiral anomaly cancellation in the lattice 
regularization (the species doubling). In the present study, we use the lattice regu- 
larization and also the Haar measure of U(N) for the average over the random vector 
potential. Therefore topologically nontrivial configurations of the vector potential 
are all included in the average. Implication of the present results for the system of 
Dirac fermions in a random vector potential with noncompact Gaussian distribution 
is discussed. 
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1 Introduction 



Random disordered systems are one of the most important problems in condensed 
matter physics. Especially in one and two dimensions effects of random disorders 
are so strong that almost all states are localized by random potentials. Nonpertur- 
bative methods are required in order to investigate random disordered systems in 
low dimensions. Recently there appeared important studies on Dirac fermions with 
random- varying mass and/or in a random vector potential. 

Study of the Dirac fermions in a random vector potential was revived in Ref.|J 
in the context of quantum Hall plateau transition. Non-Abelian generalization was 
introduced in Ref.0 in the context of a d-wave superconductor. After that, there 
appeared a number of interesting papers^, §, |, 0, §]• 

There are two technical problems for studying the system; 

(i) Normalization of states before the ensemble average over random variables. 

(ii) Integration over the white-noise random vector potential. 

For the first one, replica trick and supersymmetric(SUSY) methods are often used. 
The second problem is how to regularize the integral over the vector potential which 
has no spatial correlations. Sometimes specific parameterization of vector potential 
is used especially for the non-Abelian case|| ||. There only topologically trivial 
configurations are integrated over. By using those "technologies" , critical lines in the 
system are observed. 

Among various interesting properties of the random disordered systems, one of 
the most important problems is a disorder-induced phase transition. For the Dirac 
fermions in an Abelian random vector potential, a weak-strong disorder transition was 
found in Ref. [|J by studying multifractal scaling exponents of the critical wave function 
which is obtained exactly (see also Ref. M). Instability of the critical line is also seen 
by the existence of an infinite set of operators with negative scaling dimensions jj], 
[(J and very recently solution to the negative dimension problem was suggested 
by GurarieQ. For the case of non-Abelian random vector potential, existence of 
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operators with negative scaling dimensions was discussed in Ref . || and a termination 
mechanism was given in Ref. f| for the SU(2) case. As a related system with the 
Dirac fermion in a random vector potential, random XY model in two dimensions 
was studied in Refs. [p"0|, P| . 

In all the above previous studies, integral region of the vector potential is 
noncompact, i.e., G (—00, +00), and the probability distribution for it is taken 
to be Gaussian. However in the network models by Chalker and Coddingtonfl"2|| 
and their field-theory representation [13], the random vector potential corresponds to 
the random Aharonov-Bohm phase for electron which moves in a random potential. 
Therefore the random variables are e* Afl G U(l), and the range of the vector potential 
is compact, i.e., A^ G [— 7T, +7r]. In the weak-disorder case, the above compact distri- 
bution of the vector potential might be approximated by the noncompact Gaussian 
distribution, but at the strong disorder substantial differences in the properties of the 
system are expected to appear. 

In this paper, we shall study random Dirac fermions by employing a lattice regular- 
ization to define the system without any ambiguity. The two-dimensional(2D) Dirac 
fermions in a random vector field is formulated on the lattice and the integral measure 
of the vector potential is compact as in the original network models []12|, [L3] . There- 
fore topologically nontrivial configurations of the vector potential are all included. 
We consider U(N) vector potential for large N. Because of the compactness of the 
group U(N), the one-link integral is evaluated exactly for large N in a closed form. 
Concerning with the above problem (i), we shall use the lattice SUSY(LSUSY) meth- 
ods which we introduced before for the investigation on the SUSY gauge theory |14 



Formulation of the LSUSY is an important but still unsolved problem. However the 



LSUSY given in Ref. ||14| is suitable for the present study. 

For the one-link integral over U(N) for lage N, it is known that there are two 
regimes or "phases" , the weak-coupling and strong-coupling regimes. More precisely 
in the present context, the one-link integral exhibits a third-order phase transition as 
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the disorder strength is increased. Then one can expect that there is a genuine phase 
transition corresponding to these two regimes in the thermodynamic limit. In this 
paper we study the system in both the weak and strong-disorder cases and investigate 
the properties of both "phases" . Especially as the strong-disorder limit acquires lots 
of interest recently, we study the strong-disorder phase in detail. 

This paper is organized as follows. In Sect. 2, the model and the LSUSY are ex- 
plained. In Sect. 3, the weak-disorder regime is studied and an effective action is ob- 
tained by integrating over the random vector potential. In Sect. 4, the strong-disorder 
regime is considered. Properties of the ground state and low-energy excitations are 
clarified. It is shown that there exists a quasi-long-range order and correlation of 
the local density of states decays algebraically. Section 5 is devoted for discussion. 
Implication of the result for other interesting cases is discussed. Especially relation 
between the properties of the present model and results of the previous studies is 
examined. Physical picture of the result is explained. 



We shall study 2D Dirac fermions in a U(N) random vector potential by employing 
the lattice regularization. Action of the Dirac fermion ip a (a = 1, N) on the square 
lattice is given by 



where ) denotes lattice site, \x = (0, 1) is the direction index, U^(x) is 

U(N) field on the link (x,x + /i) U^(x) = (u^x)^ G U(N) and we set the lattice 
spacing ai = 1. The random U(N) vector potential A^ a (x) (a is the U(N) index) is 
related with U^(x) as follows, 



2 Model and LSUSY 




2 H [^l J ( x )lfiU fl (x)ilj(x + ip(x + fi)^Ul(x)if;(x)\ , 



(2.1) 



U,(x) 



(2.2) 
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where T a 's are generators of the U(N) = £7(1) x SU(N) Lie algebra. In the (naive) 
continuum limit cll — > 0, we can expand U^x) as L^x) = 1 + WLEa^^aW + 
• ■ • and recover the usual action of the Dirac fermion in the continuum. The two- 
dimensional 7- matrices are explicitly given by the Pauli matrices as 70 = <J X , 71 = a y 
and 75 = a z . 

In order to take the ensemble average over the vector potential as random vari- 
ables, we shall introduce boson field in a SUSY manner. To this end we shall slightly 
rewrite Sd in Eq.( [2.1|) . By the following transformation, 

ip(x) = T(x)x(x), ${x) = x(x)T ] (x), (2.3) 

with T(x) = ('jo) xo (li) X1 and using the identities like (7 M ) 2 = 1 and (7o)"7i = 
(— ) n 7i(7o) n {n is an integer), 

= J2x(x)D X (x), (2-4) 
where f]o(x) = 1, r]\{x) = (—) Xo , and 

Dx{x) = \Y. [ln( x ) u n( x )x(x + h)~ Vh(x - n)Ul{x - lAx{x - //)] . (2.5) 

The fields x an d X are two-component spinors but their spinor indices are diagonal 
in the action Q2.4Q . We add the following "mass term" to the action which measures 
deviation from the band center or critical line, 

S M = Mj2x(xh<jX(x) 

= MJ2[X-X+-X+X-], (2-6) 

where x = {x+iX-Y an d X = {X-jX+Y- The specific form of the above mass term 
comes from merely a technical reason which becomes clear shortly. We are interested 
in the limit M — > 0. 

We introduce a complex scalar field <fi{x) whose action is given by 

= £ D<p\x)b<p(x) + m 2 <p\x)<p(x). (2.7) 



As the modified "Dirac" operator D does not contain the 7-matrices, the same D can 
be applied for the scalar field 0. This is an essential point in the present construction 
of the SUSY lattice model. It can be shown that in the classical continuum limit the 
action in Eq.( j2.7|) reduces to the usual action of the scalar field in the continuum, 
and the integral over the scalar fields is well-defined because the action is positive- 
definite. 

The total action of the system is given by 

S = + Sff), 
S x = So + Sm 

= J2[x+DX- + X-Dx+\+Mj2[x-X+-X+X-]- (2.8) 

It is seen that the partition function is just unity if M = m for an arbitrary fixed 
configuration of the vector potential because of the cancellation of the fermion and 
boson determinants. Actually 



j [D^Dcj)} e~ s * = det-\-D 2 + m 2 

( 

J[DxD X ]e- s * = Det 



D + M 
D-M 
det(Z) 2 - M 2 ), (2.9) 



V 



where Det is the determinant of the spinor and real spaces and det is that of the real 
space. 

Moreover the action S is invariant under the following LSUSY transformation for 

M = m, 

H ] = X-e+ + X+e-, 
5 X ± = ±M<j ) e ± -D<f ) e ±} 

5 X ± = ±M0 f e± - btfe±, (2.10) 
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where e± are anticommuting spinor variables with chirality 75 e± = ±e±. 

The bosonic part of the action can be rewritten into more symmetric form 

where and are complex boson fields. By integrating over ui(x) (or <fi(x)), one can 
easily verify the equivalence of and S uv . 



The lattice Dirac action ( |2.1| ) has exact chiral symmetry for M — 0. This means 
that there appear multi-flavour Dirac fermions with opposite chirality in the contin- 
uum limit, i.e., the species doubling. S x with M = is therefore invariant under the 
following chiral U{1) x £7(1) symmetry on the lattice, 

X+(x) -> K(z)X+(X>> X-(x) -> X-(x)K! e(3;) , 

X-(x) - W e(a! )X_(x), - X+W^.), (2.12) 

where e(x) = (— ) xo+a:i and V±,W± € ?7(1). This symmetry plays a very important 
role in the discussion on the phase structure as we shall see later on. 

Expectation value of physical quantity X is given by the following functional 
integral, 

(X) = J [DUDxDxD^DtP] P[U] e~ s X, (2.13) 
where the probability distribution for the random vector potential is given by 

P[U]=exp(-Y l Tr{U li (x) + Ul{x)j), (2.14) 
9 x,fi 

[DU] = Hunk dU^x) is the Haar measure of U(N) and g is a parameter which con- 
trols the disorder strength. As we explained in the introduction, the random vector 
potentials are taken to be compact random variables. This is in contrast with the 
previous studies where the vector potential is taken to be noncompact and the prob- 
ability distribution is taken to be Gaussian. Since in the weak-disorder case, i.e., the 



case of small g , the compact distribution (|2.14j ) can be approximated as (here we 
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show the Abelian case for notational simplicity) 

e ^ cos ^, e [-7r,7r] -> e-f^", Ae(-oo,+oo) (2.15) 



5 ^0 

one may expect that the both compact and noncompact systems give similar results 
at least qualitatively. 1 However in the strong-disorder case, substantial difference 
will appear. Because of the regularization by the lattice, topologically nontrivial 
configurations are all included in the integral. This is in contrast with the discussion 
in terms of the conformal field theory(CFT)|| ||. 

In the subsequent sections we shall perform U(N) -integral in ( p. 13 ). It is known 



that there are two "phases" for this one-link integral, i.e., weak-coupling "phase" 
for small g and strong-coupling "phase" for large g. Result of the integral exhibits 
a third-order phase transition at a certain critical value of g\U§. This is merely a 
matter of kinematics of the integral over the group U(N) for large N. However we 
think that there exists a genuine phase transition in the system of the random Dirac 
fermions from weak to strong-disorder phases, i.e., disorder-induced phase transition. 



3 Weak- disorder regime 



Expectation values of physical quantities are given by ( |2.13| ) and (|2.14 ). We shall 
first perform the functional integral of the vector potential U^. Then let us consider 
the following one-link integral, 



,W(D,D) 



dU^exp[Tr(D^ + UlD, 



(3.1) 



In the present case, 



D,(x) a b 



1 Actually this expectation is too naive. The compactness of the vector potential plays a very 
important role for the correlation functions in which singular configurations of give dominant 
contribution. See discussion in Sect. 5. 
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N 



2 

Let us introduce a parameter s by 



H^Mx + V)x a (x) + <4(x + »)ip a (x) + <p\{x + AiK(x)], 
-H^-MxWix + fi) + ut(x)<p a (x + fi) + <pl(x)u a (x + fj,)]. (3.2) 



1 N 

±Y,*a 1/2 = Tr(DD)-^, (3.3) 

iV o=l 



where the x a 's are eigenvalues of j^DD. In Ref.jnj, it is shown that there are two 
regimes for the above integral (|3.1| ), i.e., weak-coupling regime for s < 2 and strong- 
coupling regime for s > 2. For D a h = ^5%, we can estimate the critical value of g as 
g c = 2 from (|3.3|). In the present case there are extra factors A and A in D and D 



and then precise value of g c cannot be determined. However we can expect that for 
sufficiently small(large) g the system is in the weak(strong)-coupling regime. 

Let us consider the weak-disorder phase first, i.e., the phase of small g. In this 



case the result of the [/(iV)-integral is given aspq| 



W(D, D) = N{2 £ atf* log^ + x\ /2 ) }. (3.4) 

a a,b 

From (U) 

(DMx)) a b = (j) 2 5 a b + jA^ b + jA a ^ + A%A^. (3.5) 

R 
N ■ 



It is straightforward to obtain (DD) 1 ^ 2 in powers of % 



—A^ + Oiig/N) 2 ). (3.6) 



Evaluation of the second term of the formula ( |3.4j ) is also not so difficult. As matrices 
D and D have both large diagonal matrix elements N/g, we can use the Taylor 
expansion for a regular function f(x,y), 

J2f( X a,Xb) = Y,f(\ + W ^-5 +W b) 
a,b a,b 9 9 



J^n\m)J { ' ^g 2, g^ Tr ^N 2 g 2 ^ N 2 g 2 ^ ' 
where io a 's are eigenvalues of (-^r — pO- Leading-order terms are obtained as, 

£ lo g(*a /2 + xl /2 ) = i(A% + Ay + \^K* A \« + 0((^/iV) 2 ). (3.7) 

a,b ~~ 

From ( |3.(j| ) and (|3.7|) , the effective theory which appears after the integration over 
the U(N) field for small g is a SUSY extension of the Gross-Neveu model.| Detailed 



study of this model will be given elsewhere ||16||. However here we mention that the 
four-Fermi coupling in the interaction term A^A^ has the sign which indicates in- 
stability of the Dirac fermion at the band center or on criticality. By the usual 1/N 
expansion, it is expected that the chiral condensation (ijjip) 7^ occurs in the present 
system. However the system under study is invariant under the chiral U(l) x U(l) 
transformation Q2.12|) , and therefore it is expected that only a quasi- long-range order 
exists as in the Kosterlitz-Thouless phase, though careful study is required because 
of the existence of the bosons. If this expectation is correct, then excitations are 
massless "pions", their SUSY partners and the Dirac fermions [ |16|1 . 

Closely related model appears in the case of 2D Dirac fermion in a random non- 
compact Abelian vector potential, i.e., a SUSY Thirring model. This is not surprising 
because in the weak-disorder limit the compact measure of the vector potential is 
well approximated by the noncompact Gaussian distribution, as we explained before. 
In Ref. |0, it is shown that there exist an infinite number of relevant operators with 
negative scaling dimensions. Very recently Gurarie suggested a possible solution to 
this problem 0. In this argument, some ad hoc cutoff regularization is used for the 
functional integral over the noncompact vector potential. On the other hand in the 
present system, the compact Haar measure is used for the integration over the vector 
potential and therefore it is expected that a cutoff appears in a natural way. We shall 
2 More precisely, Dirac fermions in this theory has "flavour" degrees of freedom as a result of the 
species doubling. 
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show that this is the case. Relation between Gurarie's argument and the present 
lattice model will be discussed later on. 

Physical picture of the above phenomenon will be discussed rather in detail in 
Sect. 5 after investigation of the strong-coupling regime. In the following section, we 
shall study the strong-disorder phase which is the main subject of the present paper. 



4 Strong-disorder regime 



In the strong-coupling regime of the £/(iV)-integral, W(D, D) is given by the following 
formula [15], 



W(D,D) 



,1/2 



1 4 

-^EM(c+^ 1/2 + (c + ^) 1/2 )}, 



a.b 



where x a 's are again eigenvalues of -^DD and a constant c is given by 



N 2 
1 

2iV 



-1/2 



(4.1) 



(4.2) 



As we explained above, the formula fl4.1|) is suitable for large g. The limit g — ► +oo 
is nothing but the strong-coupling limit of the SUSY lattice gauge theory which was 
studied in Ref.fl4|]. There we showed that the condensations like (xx)i (^V)> Re- 
occur, whereas (ip'oo) = (ou'tp) = 0. Here we assume a similar pattern of condensations 
for large g. Then it is not so difficult to calculate the effective action from ( ]4.1| ) by the 
1/g-expansion. Here again the Taylor expansion is useful to convert the summation 
over the eigenvalues x a into the trace of the matrix DD. 



After some calculation, we obtain flTT| 



—W(D,D) 



±) - F(0 - F(() 



X a {x) X a{x + ^)G{\) 



-^{x)u\{x + vl)G(£) - u\x)^ a {x + /x)G(C) 



' gN 2 



X a {x + n)xa(x)G(\) - V a (x + »H(x)G(0 



-u a (x + ^)ip{(x)G(0}+O(l/(g 2 N)) 



(4.3) 
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where A± etc are composite fields of x± etc; an d they are explicitly given by 
A± = A At± (x) = m±(x)m±(x + /i), m±(x) = -j= E X±(^)Xa T (^), 

a 

£ = = a{x)p{x + fi), C = CmO*0 = P(x)a(x + n), 

a 

Functions -F(x) and G(x) are given by 



«(*) = 4e^)^(^), /%) = 4e^)^)- ( 4 - 4 ) 



= l_(l_ x) l/2 + log[ 1 _ (1 + (1 _ x) l/2 )]) 

G(x) = (l + Cl-x) 1 / 2 )" 1 . (4.5) 

Actually there are additional terms of composites like x a { x ) l pi,{ x ), but they do not 
have nonvanishing expectation values and give only higher-order corrections in 1/N 
to the effective action of m±(x) etc. 



We expect that the 1 /^-expansion in ( [4. 3D has a finite convergence radius. Then 



it is easily verified that the effective action can be written in terms of the composites 
m(x), a(x) and /3(x), or more precisely A M (x), and Cm 0*0 -0 Then we can introduce 
elementary fields corresponding to the composite fields in the path-integral formalism. 
First for the composite "meson" field m±(x), we have identity like (up to an irrelevant 
constant), 

Z!{J) = Jd X d X e Jm = J N = J Q 2 ^(pe^- N eMJpe ie ) 

= J dM M' N exp(JM). (4.6) 

Equation ( ^4.6| ) means that the path integral of the elementary meson fields Ai± is 
defined by the above contour integral and the radius p should be taken for the angle 
integral to be well-defined, i.e., p should be a maximum or saddle point of the effective 
potential. 



'Strictly speaking, here we assume that the U(N) symmetry is not spontaneously broken. 
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On the other hand for the boson-composite field a(x), we can prove the following 
identity, 

Z^(J) = f dyd V e- Ja = J~ N = / + °°c/(ln$) $ Jv exp(-J$). (4.7) 

In a similar way we introduce elementary field ty(x) for the composite field /3(x). 

From (|4.1|) , ( |4.3|) and (|4.7|) , the effective action in the strong-disorder phase is 
obtained as 



7fS eff = - £ [F(X^(x))-\logX^(x)] -M^{M+-MJ) 

+ E [nu*)) - 7 lo § W*)] + ^E« 

^ X 

+ E [>(&(*)) - 7 log + ™ E + 0(1//), (4.8) 

x,n x 

where A M± (x) = M±(x)M.±(x + [/,), Cni x ) = &{x)^{x + and Cm( x ) = ^ / (^) < 1 ) (^ + /•*)• 
Terms of 0(1/ g 2 ) have a similar structure to the leading-order terms. 

Then it is straightforward to study the structure of the ground state and low- 
energy excitations. For vanishing masses M = m = 0, the ground state is parameter- 
ized as follows, 

, vUn±, at even sites 
(M ± (x)) = { (4.9) 

vUq±, at odd sites 

, je 01 , at even sites 
<$(*)> = { (4.10) 

ve CT2 , at odd sites 

, ve a2 , at even sites 
(*(*)) = < (4.11) 
t>e at odd sites 

where Uq± G U(l), <7j's (i = 1,2) are real numbers and v is obtained from the 
stationary condition of the effective potential, 

dV(v 2 ) dF{v 2 ) 1 



dv 2 dv 2 4t> 



2 



0, (4.12) 
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with the following solution 

v 2 = 3/4. (4.13) 

The mass terms lift the above degeneracy and determine the expectation values 
as (M.+) = —{JA-) and o\ = Oi — 0. Obviously the ground state preserves the 
SUSYQ. 

For vanishing masses, there exist the degeneracies of the ground state parameter- 
ized by Uq± and cr, which originate from the chiral symmetry fl2.12|) and its SUSY 



counterpart for (p and uj. By the Coleman- Mermin- Wagner theorem, in two dimen- 
sions continous symmetry is not spontaneously broken and there exists no long-range 
order. Therefore we cannot expect the condensations (M.±) ^ 0, etc. Instead the 
ground state exists in the Kosterlitz-Thouless phase with gapless excitations. In fact 
we can explicitly show the existence of massless modes which destroy the off-diagonal 
long-range order. These excitations are described by the "pion" fields, 

!vU±(x) = ve ln± ( x \ at even sites 
(4-14) 
vU±(x) = ve ,7r ±w, at odd sites 

and similar SUSY excitations for and *f?(x). From ( |4.10|) and (|4.11|) , it is obvious 
that these low-energy excitations are nothing but "density wave" of the SUSY bosons 
which is commensurate with the lattice structure. 

Effective action of ir±(x) is obtained as follows from (|4.8|) and (|4.14 ), 



S w = yC^[7r ± (a; + /i) -n±(x)] 2 , 

C = F"(v)v 4 + i, v 2 = 3/A. (4.15) 

Therefore the correlator of JA±(x) exhibits a power-law decay 

(^(x)^(0)) = (M±(x)M±(0)) ~ \x\- 1/{2wNC) . (4.16) 

Corrections of 0(1/ g 2 ) can be calculated systematically and the scaling dimension of 
A4± acquires correction of 0{l/g 2 ). 
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There appear no signs of instability of the ground state. Low-energy excitations 
in the boson sector are given by local fluctuations of <j{ (i = 1,2) in Eqs.( j4.10| ) and 
( |4.11| ). They are SUSY counterparts of ir±(x) and stable. All correlation functions 
in the fermion sector have nonsingular behaviour like Eq.( 4.1C ). This is in sharp 
contrast with the previous results which show that fermion composite operators (as 
well as boson composite operators) have negative scaling dimensions |], H ||. We 
think that the compactness of the functional-integral measure of the vector potential 
plays a very important role for the stability. This important point will be discussed 
in the following section. 

From the discussion given so far, it is obvious that the algebraic decay of the cor- 
relation functions in the strong-disorder phase comes from the exact chiral symmetry 
on the lattice which is a result of the species doubling. In the single-flavour case in 
the continuum, there exists anomaly in the chiral symmetry because of the coupling 
with the vector potential, and therefore the genuine condensation (ijjip) 7^ is possible 
even in two dimensions just as in the Schwinger model.0 



5 Discussion 

In this paper we studied Dirac fermions in a random U (N) vector potential. We em- 
ployed the lattice regularization and the compact Haar measure in order to make the 
functional integral over U (N) vector potential well-defined. In this formalism topo- 
logically nontrivial configurations are all integrated over. This is in sharp contrast 
with the approaches given so far. The ensemble average over the random vector po- 
tential was taken by introducing bosons in a SUSY way. We think that this approach 
is important because there appeared some evidences that there exists a disorder- 
induced phase transition in the present system. In order to investigate this problem, 
In other words, the would-be Nambu-Goldstone boson acquires a mass by the chiral anomaly, 
and it generates no severe infrared singularities. 



15 



a well-defined formalism is indispensable. 

For the one-link U(N) integral, it is known that there are two regimes, i.e., the 
weak and strong-coupling regimes, which correspond to the weak and strong-disorder 
cases, respectively. We obtained effective theory by integrating over the vector po- 
tential in both regimes. In the weak-disorder phase, the effective theory is a SUSY 
extension of the Gross-Neveu model. We call this regime phase A. Detailed studies 
on the effective field theory of the phase A will be reported elsewhere, but sign of the 
effective coupling constant of the four-Fermi interaction indicates instability of the 
ground state to the state with the chiral condensation. From the investigation of the 
Gross-Neveu model by the 1/N expansion, we expect that there appears the chiral 
condensation with a quasi-long-range order. 

Then we studied the strong-disorder phase rather in detail. We call this regime 
phase B. We showed that in the phase B the density operator of the fermion has the 
quasi-long-range order and low-energy excitations are the "pions" and density wave 
of the bosons whereas no Dirac fermions with the original U(N) quantum number 
appear. This result stems from not only the strong-disorder properties of the vector 
potential but also the exact chiral symmetry and its SUSY counterpart on the lattice, 
i.e., anomaly cancellation by the species doubling. Therefore, the result indicates that 
genuine condensation of the fermion density operator occurs in the single-flavour case 
with chiral anomaly just as in the Schwinger model. Anyway, we expect the existence 
of the phase transition at a certain critical value of the disorder strength g c from the 
phase A to B.Q 

Let us discuss physical picture of the above phenomena. Coupling with the vector 
potential reduces the effective hopping of fermions which is simply given by t e g = 
t ■ (U^x)) where t is the original hopping parameter. Obviously as increasing the 
disorder strength g, t e Q decreases and fermionic states tend to localize. Study in 
this paper shows that in the present system that localization phenomenon occurs 
5 Order of this phase transition can be of third-order as the one-link integral indicates. 
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with the chiral condensation of fermions which makes the fermions massive. As g 
is increased further more, even a local hopping of a single fermion cannot occur 
anymore and movement of a fermion always accompanies the same movement of an 
anti-fermion for fluctuations of the vector potential in the hopping cancels out with 
each other, i.e., (U^(x)) = but (U tl (x)W t (x)) = 1. That is, only a "bound state" of 
fermion and anti-fermion pair can move in this phase. Then condensation of fermion 
density operator is generated. This phase (the phase B) is more or less similar to 
the conventional confinement phase of the strong-coupling gauge theory. From this 
picture, disorder-induced phase transition is naturally understood. 

Finally let us discuss relation between the results in this paper and previous stud- 
ies. For both the Abelian and non-Abelian cases, it is known that there exist an 
infinite number of relevant operators with negative dimensions if the noncompact 
Gaussian distribution is used for the vector potential. It indicates some instability 
of the critical line. On the other hand in the present study, no signs of instability 
appear at least in the effective action obtained by integrating over the vector poten- 
tial!] We think that this stability stems from the compactness of the integral measure 
of the vector potential. Actually as recently Gurarie discussed^] and the study on 
the random XY model shows |TlJ|, the instability comes from the noncompactness of 
the vector potential. More explicitly in the discussion in the continuum, the vec- 
tor potential is parameterized as follows (we here consider the Abelian case for 
simplicity^) , 

Au = ^ud v 6 + <9 M ?7, e i = -eio, (5.1) 

where 6(x) and r)(x) are scalar fields and 8(x), r](x) G (— oo, +oo). The instability and 

the negative dimensions of the relevant operators essentially come from the following 
6 Strictly speaking, the correlation function in the boson sector (^(x)^f(0)) tends to diverge for 

\x\ — > oo. However correlators in the fermion sector, which are physical quantities in the present 

system, exhibit no singular behaviour. 

7 For the Abelian case, similar discussion of the strong-disorder case is possible. See Ref . Q . 
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correlation function of 0(x)||ll||, 



where c is a real number and the above expectation value is evaluated with the 
following probability distribution, 

P[9] oc exp { - - f d 2 x {d„9) 2 }. (5.3) 

Then it is not difficult to show that the operator e cd ^ has a negative dimension and 
the correlator in Eq. (|5.2|) tends to diverge for large \x\ — x%\. However it is obvious 
that if the integral region of 9{x) is compact, this divergence does not occur. This is 
an essential point of Gurarie's argument || . 

Gurarie used some ad hoc cutoff regularization for the functional space of 9(x). 
On the other hand in this paper, we use the Haar measure which is compact, and 
therefore we expect that a similar cut off appears naturally in the present formalism. 
Actually we can parameterize the vector potentials U^x) in terms of two U(l) fields 
u(x) and v (x + ^±i) (6(1) is the unit vector of the 0(1) direction), where u(x) is 
defined on the sites of the original lattice and v(x + 2±1) j s on the sites of the dual 
lattice, 

U Q (x) = u(x + 6)u*(x)v(x + ~^—)v*(x + — g— ), (5-4) 

and similarly for U\(x) where u*(x) is the complex conjugate to u(x) etc. It is not 
difficult to show that if we impose the conditions like Y\ x u = Yl x v — 1, then there is 
no ambiguity in this parameterization. It is obvious that 8(x)(r](x)) in the continuum 
is related to v (x + ^^)(u(x)) in ( |5.4j) as follows, 

6 ~ lnf, i] ~ ln-u. (5-5) 

Therefore the integral region of 9(x) is compact in the present formalism, 9(x) G 
[— 7r, +7r], i.e., there exists the natural cutoff. Please notice that this compact region 
of 9{x) remains the same even if we recover the lattice spacing ol. 
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expression (5.1 



In the sense explained above, the model in the present paper is different from 
those with noncompact random vector potential which were studied in the previous 



papers. Our model is close to the network models by Chalker and Coddington|T2] 



and their field-theory models [13]. Then it is not so surprising that the stable ground 
state appears and the low-energy excitations are the massless "pions" etc as in the 
Kosterlitz-Thouless phase even at the strong-disorder limit. Also we can conclude 
that the disorder-induced phase transition which we found in this paper is a new one. 
This phase transition is expected to be of a topological nature of the vector potential 
but it is not definitive at this stage. 
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Note added 

After submitting this paper, we got acquainted with the paper by Altland and 
Simons [|H^[ which also studies the random flux model on a lattice. 
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